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ABSTRACT 

We report the results of a sensitive search for 12. 2- GHz raethanol maser emission 
towards a sample of eight high-mass star formation regions in the Large Magellanic 
Clouds which have been detected in other maser transitions. We detected one source 
towards the star formation region N105a. This is the first detection of a 12.2-GHz 
methanol maser outside our Galaxy. We also made near-contemporaneous observations 
of the 6. 7- GHz methanol and 2 2- GHz water masers towards these sources, resulting 
in the detection of water maser emission in six new sources, including one associated 
with the strongest 6.7-GHz maser in the MageUanic Clouds /i?^S'05011-6815. The 
majority of the maser sources are closely associated with objects identified as likely 
Young SteUar Objects (YSO) on the basis of Spitzer Space Telescope observations. 
We find that the YSOs associated with masers tend to be more luminous and have 
redder infrared colours than the sample as a whole. SED modeling of the YSOs shows 
that the masers are associated with sources of higher central mass, total luminosity 
and ambient density than the majority of YSOs in the LMC. This is consistent with 
the well-established relationship between luminous methanol and water masers and 
young, high-mass objects observed in the Galaxy. 
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1 INTRODUCTION 

The Magellanic Clouds represent a unique laboratory for 
studying star formation in an environment which differs 
from that of our Gala:xy. They are nearby (in comparison 
to most other galaxies), at distances of approximately 50 
and 60 kpc respectively for the Large and Small Magel- 
lanic Clouds (Feast 1999; Walker 1999). This enables them 
to be studied with intrinsically high sensitivity and resolu- 
tion. The relative star formation rates in both the Large 
Magellanic Cloud (LMC) and the Small Magellanic Cloud 
(SMC) are greater than in the Milk y Way, but th e metal- 
licity is significantly lower (jRussell fc Dopitalll992i l. These 
conditions are thought to be similar to those in larger galax- 
ies earlier in the history of the Universe. Hence studies of 
star formation in the Magellanic Clouds can yield insights 
into star formation at redshifts of jz > 1. 
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Interstellar masers are excellent signposts of high-mass 
star formation in the Milky Way. They are intense and oc- 
cur at wavelengths that are not si gnificantly absorb ed by 
intervening clouds of gas and dust jcreen et al.ll2009l '). The 
first complete survey of the LMC and SMC was recently 
undertaken in the 6.7-GHz transition of meth anol and the 
6.035-GHz transition of OH (jCreen et al.ll2008l '). The survey, 
made with the methanol multibeam receiver on the Parkes 
telescope, had a 1-a sensitivity of 0.13 Jy and 0.09 Jy over 
the entire SMC and LMC respectively (the extent of both 
Clouds, as defined by the CO and/or Hi emission). The re- 
gion of the LMC containing the bulk of the CO emission and 
all the previously known masers were surveyed to a better 
1-cr sensitivity of approximately 60 mjy. Despite this com- 
prehensive and sensitive search only four 6.7- GHz methanol 
masers have been discovered in the LMC l|Sinclair et al 



19921 : lEllingsen et"ai]|l994l : iBeaslev et al]|l99a I Green et al 



20081). Taking into account the relative rates of high- mass 
star formation in the two systems, this is a factor of 4-5 lower 
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than expected on the basis of the observed lu minosity dis- 
tribu tion of Galactic 6.7-GHz methanol masers (| Green et alj 
l2008h . In contrast, the number of detected ground-state OH 
and water masers is roughly consistent with the difference in 
the star formation rates between the Milky Way and LMC. 
This discrepancy has been attributed to the lower metallic- 
ity of the LMC. 

The current rate of high-mass star formation within 
the Magellanic clouds, and its relationship to the dynamical 
interaction of the clouds with the Milky Way, are a topic 
of active investigation. The Spitzer Space Telescope Legacy 
program SAGE (Surveying the Agents of a Galaxy's Evolu- 
tion) has made observations of a number of nearby galaxies 
(including the LMC) to exam ine the current star forma- 
tion rate (jMeixner et al.l [20061 ). The unique association of 
6.7- and 12.2-GHz methanol masers with young, high-mass 
star forming regions has made them very useful signposts 
of such regions within our Galaxy. A recent comparative 
study of 6.7- and 12.2-GHz methanol masers in the Milky 
Way has shown that 12.2-GHz masers are pref erentially as- 
socia te d with more luminous 6.7-GHz sources (|Breen et al.l 
l2009l ). [ Breen et al.l also used their observations and com- 
parisons with published 1.2-mm dust continuum and mid- 
infrared data to st art to quantify the m aser evolutionary 
scheme proposed bv lEUingsen et al] (|2007l ). They found that 
the 12.2-GHz methanol masers are associated with a later 
evolutionary phase than those sources which show emission 
from only the 6.7-GHz transition. The median ratio of the 
peak flux density of the 12.2- to the 6.7-GHz transitions 
was found to be 1:5.9, similar to previous studies. It was 
also shown that the ratio depends on luminosity (and hence 
age), increasing as the luminosity o f the 6.7-GHz emission 
increases (see Fig. 4 of lBreen et al.h 

In order to investigate whether the properties of 12.2- 
GHz methanol masers in the LMC differ from their Galactic 
counterparts we have undertaken new, sensitive observations 
towards kno wn 6.7-GHz methanol and OH maser sources 
in the LMC. lEUingsen e\~^ l|l994) searched for 12.2-GHz 
methanol masers towards the two sources known at the time 
Nil (MC18) and N105a (MC23), with 3-cr detection upper 
limits of 0.1 and 0.3 Jy respectively. There are no published 
observations at 12.2-GHz towards the two more recently dis- 
covered 6.7-GHz methanol masers in the LMC IRASQhQll- 
6815 and N160a (MC76), nor at the locations of any of the 
known OH or water masers. 

In contrast to methanol masers, the abundance of 22- 
GHz water masers in the Magellanic clouds appears to be 
consistent with that observed in the Milky Way. This sug- 
gests that water masers may be more useful signposts of 
young high-mass star formation regions within the LMC. 
In order to better characterize the properties of methanol 
and water masers in the LMC we have also made observa- 
tions of the 6.7-GHz methanol and 22-GHz water masers 
in the LMC which were nearly contemporaneous with our 
12.2-GHz search. 



2 OBSERVATIONS & DATA REDUCTION 

We have searched a sample of eight known maser regions 
in the LMC for the presence of 12.2-GHz methanol masers. 
The sample included all known ground and excited-state OH 



masers, 6.7-GHz methanol masers and 22-GHz masers in the 
Magellanic Clouds. Observations were carried out with the 
Australia Telescope National Facility (ATNF) Parkes 64-m 
radio telescope during 2008 June 20-25. The observations 
were made with the Ku-band receiver which detected two 
orthogonal linear polarizations. It had typical system tem- 
peratures of 205 and 225 Jy for the respective polarizations 
throughout the observations. The Parkes multibeam correla- 
tor was configured to record 8192 spectral channels across a 
16-MHz bandwidth for each polarization. This corresponds 
to a velocity coverage of ~290 km s~^ and after Hanning 
smoothing, a spectral resolution of 0.08 km s~^. The Parkes 
radio telescope has RMS pointing errors of ~10 arcsec and 
at 12.2-GHz the telescope has a half-power beam-width of 
1.9 arcmin. 

All sources were observed at a fixed frequency of 12167- 
MHz (i.e. with no Doppler tracking), which alleviated the 
requirement for a unique reference observation to be made 
for each of the source positions. The data were processe d 
using the same method as described in iBreen et al.l l|2009l ). 
except that a running median of 50 channels (as opposed 
to 100 channels) was used to correct the spectrum baseline. 
The data were reduced using the ASAP (ATNF Spectral 
Analysis Package). Alignment of velocity channels was car- 
ried out during processing. Absolute flux density calibration 
was achieved by observing PKS B1934-638 each day which 
has a n assumed flux density of 1.825 Jy at 12167-MHz (|Saultl 
l2003h . Each of the fields was observed for 160-170 minutes, 
resulting in a typical RMS in the final spectrum (after Han- 
ning smoothing and averaging over polarizations and time) 
of 27 mjy and a corresponding 5-(j limit of .135 Jy. The 
adop ted rest frequency was 12.178597 GHz l|Muller et al.l 
l2004h . 

The 6.7-GHz observations reported by ICreen et al.l 

(|2008l ) were made in 2006 June and September, approxi- 
mately 2 years prior to the 12.2-GHz observations reported 
here. As masers are known to be variable we undertook ad- 
ditional 6.7-GHz observations using the Parkes 64-m radio 
telescope in 2008 August (< 2 months after the 12.2-GHz 
observations) using the methanol multibeam receiver. The 
observations were made in MX mode, with 10 minutes per- 
beam on source toward each of the f our k nown 6.7-GHz 
methanol maser sites (see lGreen et al. I l|2009t ) for a detailed 
description of the system parameters) . We recorded a 4 MHz 
bandpass for each beam, for each of two circular polariza- 
tions for both the 6.7-GHz methanol and the 6.035 MHz 
OH transitions. The spectra were produced by an autocor- 
relation spectrometer with 2048 spectral channels. The data 
were not Hanning smoothed and the resulting velocity res- 
olution for the 6.7-GHz observations was 0.11 km s-^ This 
is very similar to the velocity resolution of the 12.2-GHz ob- 
servations. A fault with one of the beams (beam 3) meant 
that its data had to be excluded, resulting in an on-source 
integration time of 60 minutes per source. 

The 22-GHz water maser observations were made with 
the Australia Telescope Compact Array (ATCA) during 
2008 August 18,19, with the array in the 6B configuration. 
The observations were made in Director's time allocations 
and spanned approximately 7 hours on the first day and 5 
hours on the second day. A total of eight different pointing 
centres were observed. These covered the locations of all the 
11 previously known star- formation maser sites in the LMC. 
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The observing band was centred at a fixed frequency of 
22213 MHz. The correlator was configured with 1024 chan- 
nels across a 16 MHz bandwidth recording a single linear 
polarization and also allowing measurement of continuum 
emission. The data were not smoothed, giving a velocity 
resolution of approximately 0.25 km s~^ and a usable veloc- 
ity range of approximately 195 km s~^. PKS B1934-638 was 
used as the primary flux calibrator, which at t his frequenc y 
is assumed to have a flux density of 0.83 Jy (jSaultl [20031 1. 
Observations of PKS B0537-441 were used to calibrate the 
bandpass and PKSB0637-752 was used as the phase cali- 
brator. The observations were made as a series of 90 second 
integrations. These were regularly interleaved with obser- 
vations of the phase calibrator. In total, each source was 
observed 9 or 10 times over the two observing sessions re- 
sulting in an onsource time of 13.5 or 15 minutes per source. 
We also used the line-free spectral channels from the water 
maser observations to search for continuum emission associ- 
ated with the star formation regions. The observations had 
only a 16-MHz bandwidth and thus were not very sensitive; 
however they do have higher spatial resolution than previous 
studies of these regions. The typical RMS for the continuum 
images was 0.2-0.3 mjy beam~^. 



3 RESULTS 

3.1 12.2-GHz methanol masers 

The known LMC masers are contained in eight distinct fields 
which were targets for a 12.2-GHz methanol maser search. 
Table [T] lists either the detected peak fiux density and veloc- 
ity, or the 5-a noise level in the final spectrum. A heliocen- 
tric velocity range of 120 - 450 km s~^ was observed for each 
source. The only source detected was N105a (MC23), which 
coincidentally was also the first 6.7-GHz methanol maser 
detected in the LMC. The 12.2 GHz spectrum of N105a is 
shown in the upper panel of Fig. [T] The peak of the 12.2- 
GHz emission is 0.12 Jy (4.6-cr). As it is within the velocity 
range of the 6.7-GHz emission, and is also apparent when 
the total 12.2-GHz dataset is split in half, or when the two 
polarizations are examined independently, we are confident 
that it is a genuine detection. 

Comparing our 12.2-GHz spectrum of N105a with the 
6.7-GHz emission observed in 2006 June ijGreen et al.ll200^ 1 
shows that the 12.2-GHz peak coincid es with the s econd 
strongest spectral peak in this source. iBreen et all (|2009l ) 
found that for approximately 75 per cent of the sources 
where they detected both 12.2- and 6.7-GHz emission the 
peak of the two spectra occurred at the same velocity. For 
the remaining 25 per cent the 12.2-GHz peak coincided 
with a secondary 6.7-GHz feature. Although the associa- 
tion of the 12.2-GHz methanol maser with a secondary 6.7- 
GHz spectral feature is not unusual, it is also possible that 
the 6.7-GHz maser may have varied in th e period between 
the observations reported by ICreen et al.l (June 2006) and 
these observations. Comparison of the August 2008 6.7-GHz 
methanol m aser spectrum, with that shown in Fig. 3 of 
I Green et al.l ([2008 ) shows no detectable variation over the 
two-year period between the observations. The observations 
of the other known 6.7-GHz methanol masers in the LMC 
and the 6.035-GHz OH maser in N160a also show no sig- 
nificant changes in the absolute, or relative intensity of the 
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Figure 1. The top plot shows the 12.2-GHz emission in N105a. 
The maser peak is at 250.8 km s~^ (Heliocentric velocity frame), 
with a peak 4.6 times than of the RMS noise in the spectrum. 
The 12.2-GHz peak coincides in velocity with a secondary peak 
in the 6.7-GHz spectrum and a comparison of the two spectra is 
shown in the bottom plot with the 12.2-GHz emission shown as 
a solid line and the 6.7-GHz as a dashed line (NOTE; a reduced 
velocity range is shown on the comparison plot to allow greater 
detail to bo seen). 



maser features compared to those published bv lCreen et al.l 
(see Figure [2]). 



3.2 Water maser observations 

The current water maser observations are summarized in Ta- 
ble [l] The previously published searches for OH, water and 
methanol masers in star formation regions in the LMC span 
a period of more than 25 years. As all the maser observations 
reported in this paper were made within a 10 week period 
between mid- June and the end of August 2008, they pro- 
vide the first near-contemporaneous observations of the ma- 
jority of star formation maser transitions in the LMC (only 
the ground-state OH masers were not observed). Figure [2] 
shows the emission from the 6.7-GHz methanol, 6.035-GHz 
OH and 22-GHz water masers detected towards the known 
star formation maser sites in the LMC. This represents all 
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Table 1. A summary of the maser observations towards 12.2- and 6.7-GHz methanol masers and 22-GHz water masers undertaken as 
part of this project. All velocities are given in the heliocentric frame of reference. For the 12 sites where water masers were detected, the 
position is that of the water maser emission observed in our AT CA observations. For the remaining four sites the position is that of the 
methanol or OH maser emission reported bviGre en et al.l | |2008| ). New detections are indicated by an asterisk after the peak flux density 
for that transition. Where observations were made, but no emission was detected the S-cr limit for the spectra are given (ATCA spectra 
were formed from the image cubes). 
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the known star formation masers in the LMC with the excep- 
tion o f the 1.6-GHz OH transitions fsee lBrooks &: WhiteoakI 
Il997l . for spectra of the OH masers). 

Sensitive, high resolution water maser observations to- 
wards known or suspected high-mas s star formation region s 
have previously been published by iLazendic et~all l|2002h . 
Their observations were also made with the ATCA, but us- 
ing only three antennas and a maximum baseline length 
of 200m. We detected all but one of the sources they ob- 
served and found six new water maser sources. The excep- 
tion was in the N113/MC24 field, where they reported the 
presence of a water maser at 05''13'"17.67^ -69°22'21.4" 
(J2000) with emission in the velocity range 247.3 - 254.4 
km s~^. We did not detect any emission from that loca- 
tion, but did detect two new water masers in this field. The 
first is at 05''13™ 17.17", -69°22'21.6" (J2000), offset by 2.6 
arcseconds, nearly entirely in right ascension from the non- 
detected site. It has a peak velocity at 257.1 km s^^ and a 
much broader velocity range than any of the other sources 
seen in the region. The second has a velocity range slightly 
lower than the non-detected source, peaking at 249.1 km s~^ 
at 05''13'"17.69", -69°22'27.0" (J2000) (i.e. it is six arcsec- 
onds further south, but with essentially the same right as- 
cension). The locations of the various masers, and infrared 
sources in the N113 region, along with the associated 22- 
GHz radio continuum emission is shown in Figure [31 We 
ha ve confirmed wi th the authors that the position reported 
Lazendic et al.l is correct (John Whiteoak, private com- 
munication) and it appears likely the non-detection is sim- 
ply due to the variability often exhibited by water masers. 
The infrared and radio continuum emission suggests that the 
newly detected water masers (and the non-detection) may be 
associated with an energetic outflow from the nearby YSO. 
The only newly detected water maser in a region where 



this species had not previously been detected was towards 
/i?AS'05011-6815 (Fig. [gj. This is the site of the strongest 
6.7-GHz methanol maser in the LMC, which has a peak 
flux density of approximately 3.5 Jy. We detect water maser 
emission with a peak flux den sity of 0.08 Jy toward the 
location of the methanol maser. iLazendic et"all (|2002l ) also 
searched for water maser emission towards this source, but 
did not detect any emission (with a limit of 0.3 Jy) at the 
epoch of their observations. 

lOliveira et all (|2006l ) used the Parkes 64-m telescope 
to search for water masers in several regions in the LMC 
and identified several new water masers in the 30 Doradus 
(N157) region. Our ATCA observations confirm the presence 
of multiple water masers in this r egion, altho ugh their rela- 
tionship to those of lOliveira et al.l is not clear. [Oliveira et al.l 
(2006) reported four masers in the N157 region, which they 
labelled 0539-691a, b, B and C. We detect three masers, but 
all are within 35 arcseconds of the stronge st maser in the 
region (0539-691a) previously observed by ILazendic et al.l 
(2002). Their source 0539-691B lies beyond the half-power 
point of the primary beam of our ATCA observations, while 
0539-691C is close to the edge of the beam and hence it is not 
surprising that we d o not detect either of these sources. The 
position reported by lOliveira et all (|2006t ) for the strongest 
maser 0539-691a is offset 21.3 arcseconds (18.1 arcsec later 
in right ascension and 11.3 arcsec further south in declina- 
tion) from the position we have determined with the ATCA 
(which is les s than an arcsec ond different from the position 
reported by ILazendic et alj ) . However, applying the same 
pointing correction to the position reported for 0539-691b 
does not improve its agreement with the location of either 
of the new maser sources we have detected. These differences 
most likely just reflect the difficulty of measuring the rela- 
tive location of masers with overlapping velocity ranges with 
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Figure 2. Maser spectra of the 6.035-GHz OH, 6.7-GHz methanol and 22-GHz water masers associated vi^ith star formation regions in 
the LMC (except for N105a/MC23 which is shown in Figure 1). Where multiple water masers were detected towards one region, the 
order of the spectra follows the order of the source positions in Table ^ The spectra are ail from the Parke s or ATCA observa tions 
described in this paper, with the exception of the 6.035-GHz OH spectrum for N157a/MC74 which is taken from ^ Green et ah Il l2008l'l . For 
the 6.035-GHz OH spectra the solid Une is right-hand circular polarization and the dashed line left-hand circular polarization, and the 
flux density scale refers to a single polarization. NOTE: The negative spike at 253 km s"^ towards the N160a/MC76 at 05''39'"43.83* 
-69°38'33.8" is due to sidelobe interference from the nearby source at 05''39™38.97'' -69°39'10.8". 
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Figure 2 - continued 
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Figure 3. The distribution of four water masers (crosses) in the N113 region. The stro ngest water maser in t he LMC is located in the 
bottom left of this image, while the asterisk marks the location of the maser detected bv lLazendic et ahl ||2002| ) which we did not detect. 
The three black circles mark the location of YSOs identified bv lGruendl fc Chul ||2009|'| . The greyscale image and contours both show the 
22 GHz radio continuum emission, the contours are at 3, 6 mjy beam~^. 
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Figure 4. Tlie distribution of water mas ers (cross) and tlie loca - 
tion of tiie YSO (black circle) identified bv lGruendl fc Chul ||2009| ) 
in the N159 region. The greyscale image and contours both show 
the 22 GHz radio continuum emission, the contours arc at 3 and 
4 mjy beam~^. 



separations less than the primary beam using a single-dish 
telescope. 



3.3 22-GHz radio continuum emission 

We detected weak 22-GHz continuum emission in three fields 
N113/MC24, N159 and N160a/MC76, with peak fluxes of 
12.9, 2.7 and 42.5 mJy beam'^ For N159 and N160a/MC76 
the radio continuum sources are coincident with the location 
of the water masers to within 0.1 arcseconds. The relative 
distributions of the water masers, infrared YSO and radio 
continuum emission are shown in Figs|3]-[5l For fields where 
no emission was detected the 5-cr limit was typically around 
3 mJy beam~^, with the highest being 5.2 mJy beam~^ for 
N105a/MC23. 



4 DISCUSSION 

4.1 Maser abundances in the LMC 

Previous studies have noted the apparent underabundance 
of methanol masers in the Magellanic C l ouds in compariso n 
to the Milky Way (|Beaslev et all 1 19961 : ICreen et al.[|2008l ). 
The current observations represent the first search for 12.2- 
GHz methanol maser towards a complete sample of 6.7-GHz 
methanol and 6.035 GHz OH masers. In addition to covering 
all the known 6.7-GHz methanol and 6.035 GHz OH masers 
in the LMC, it also covers all the known ground-state OH 
sources, and 22-GHz water masers associated with a star 
formation region. 

There are a number of publications which give a detailed 
comparison of the properties of large samples o f Galactic 6.7- 
and 1 2.2- GHz methano l mase rs, in particular ICaswell et al.l 
l| 19951 ) and lBreen et all (|2009l ). Figure [6] shows a log-log plot 



of the ratio of the 6.7- to 12.2-GHz peak versus the 6.7- 
GHz isotropic luminosity (determined from the peak flux 
density and measured in units of Jy kp c'^) for a sam ple of 
63 Galactic 12.2-GHz methanol masers (|Breen et al.l i2009). 
Figure |6] includes the four LMC points as fllled squares. 
Our upper limits for the 12.2-GHz emission towards the 
6.7-GHz methanol masers other than N105a, allows us to 
determine a lower limit to the 6.7- to 12.2-GHz peak flux 
ratio for these sources. These lower limits are indicated by 
the presence of an upward pointing arrow. Fig. [6] shows that 
the ratios (or limits) for 6.7- to 12.2-GHz emission in the 
LMC lie within the range observed in Galactic sources. The 
dashed vertical line indicates the 5-a det ection limit for th e 
methanol multibeam survey of the LMC (|Green et al-lbOOSl ). 
The dashed sloped line indicates how the lower-limit on the 
6.7- to 12.2-GHz peak flux density ratio changes as a func- 
tion of 6.7-GHz peak luminosity for sources at the distance 
of the LMC for observations with a 5-a detection limit of 
0.135 Jy (the limit for our 12.2-GHz obse rvations). Ou r ob- 
servations, when combined with those of I Green et al.l show 
that there are no methanol masers in the LMC with charac- 
teristics which place them in the lower-right region of Fig. [6] 
(as delineated by the dashed lines). However, only 20 per 
cent of Galactic masers with both 6.7- and 12.2-GHz fall 
within these limits (which denote the region occupied by lu- 
minous 6.7-GHz methanol masers with a relatively strong 
12.2-GHz counterpart). Two of the lower- limits established 
for the 6.7- to 12.2-GHz peak flux density ratio are well 
below the median value observed in Galactic sources for 
6.7-GHz masers of comparable luminosity. However, obser- 
vations approximately two orders of magnitude more sen- 
sitive would be required to establish whether these sources 
are under-luminous at 12.2-GHz compared to similar Galac- 
tic sources. For the most luminous LMC 6.7-GHz methanol 
maser /7?^S'05011-6815, the established lower limit shows 
that any 12.2-GHz methanol maser in this source is less lu- 
minous than for the majority of similar Galactic 6.7-GHz 
methanol masers, but well within the observed range. In 
summary, with our existing data there is no evidence for any 
significant difference between the ratios of methanol maser 
lines seen in the LMC compared to the Milky Way. Fig. [S] 
also nicely illustrates that the majority of the known 6.7- 
GHz methanol masers in the LMC are not particularly lu- 
minous in comparison with the luminosities of the Galactic 
sample shown here. 

We can also compare the rate of association and relative 
strength of the 12.2-GHz methanol and ground-state, main- 
line OH masers with those observed in Galactic samples. 
ICaswelj 1 199 8') provides the most comprehensive, sensitive 
sample of Galactic ground-state OH masers in star forming 
regions, 80 per cent of which have an associated 6.7-GHz 
methanol maser. There are a total of four g round-state OH 
masers known in the LMC (|Brooks fc Whit eoak 1997), and 
on the basis of Galactic statistics we would expect 3-4 of 
these to have associated 6.7-GHz methanol masers, b ut only 
two do (/7?^505011 & N160a). iBreen et all (|2009l ) found 
that 12.2-GHz methanol masers are found preferentially to- 
wards 6.7-GHz methanol masers with an associated ground- 
state OH masers. Our single detection does not follow this 
trend. 

In sources where 1665-MHz OH, 6.7-GHz methanol 
and 12.2-GHz methanol are all observed, both the OH 
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Figure 5. The distribution of two wa ter masers (crosses), an d the methanol masers (square) in the N160a region. The six black circles 
mark the location of YSO identified bv lGruendl fc Chul | |2009|) . The greyscale image and contours both show the 22 GHz radio continuum 
emission, the contours are at 9 and 26 mjy beam~^. 



and 12.2-GHz methanol are, typically , significantly weake r 
than the 6.7-GHz methanol masers (jCaswell et al.lll995h . 
ICaswell et al] l|l993h studied the relative intensity of 12.2- 
GHz methanol and ground-state OH lines in a sample of 173 
OH masers. They found a large scatter in the observed ratio, 
with a 12.2-GHz:OH median ratio of 5:1 for the 53 sources 
where both transitions were detected, falling to 1:1 if 12.2- 
GHz non-detections were included in the sample. This sug- 
gests that where both transitions are observed the 12.2-GHz 
methanol masers will usually be stronger. 

The small sample size for all the different maser tran- 
sitions observed in the LMC, along with the wide range of 
intensity ratios observed in Galactic samples does not allow 
us to draw any firm conclusions as to whether the proper- 
ties of the LMC masers differ from those seen our Galaxy. 
However, the absence of 12.2-GHz methanol masers towards 
sources with both 6.7-GHz methanol and ground-state OH 
masers, and the rate of detected 6.7-GHz methanol towards 



OH masers in the LMC are both consistent with the estab- 
lished trend for the Magellanic Clouds to have an under- 
abundance of methanol masers compared to the Galaxy. 

Comparison of the relative rates of association between 
methanol and water masers is more difficult because less 
work has been done on this in the Milky Way. At present 
there are no published large-scale unbiased water maser 
surveys to compare to the methanol maser surveys. How- 
ever, we note that this will change in the relatively near 
future when the results of the H2O galact ic Plane Survey 
(HOPS) are published (| Walsh et al.l I2OO8I ). However, un- 
til the results of HOPS are available, the only information 
we have comes from tar geted s earche s of one tra nsition to- 
wards the other ( e.g. | Szvmczak et al.il2005l : IXu et al. 200^). 
ISzvmczak et al.l (|2005[) found a detection rate of approxi- 
mately 50 per cent for water masers towards k nown 6.7-GHz 
methanol masers. Conversely IXu et al.l l|2008l ) detected 6.7- 
GHz methanol masers towards approximately 30 per cent of 



Masers in the LMC 9 




1 2 3 



log{6.7-GHz peak luminosity) 

Figure 6. A plot of the ratio of tiie 6.7- to 12.2-GHz peak 
flux density versus the peak luminosity of the 6.7-G Hz methanol 
maser. The open circles are Galactic sources from iBreen et alj 
(|2009l1 . the filled squares are the LMC measurements, with those 
which are lower limits (i.e. all but N105a), indicated by upward 
arrows. The solid line is the linear regression relationship between 
6.7-GHz peak luminos ity and 6.7- to 12.2-GHz peak ratio deter- 
mined by iBreen et al.l for Galactic sources. The dashed vertical 
line represents the 5-ct detection limit for the methanol multi- 
beam survey of the LMC (assuming a distance of 50 kpc). The 
dashed sloped line shows the lower-limit on 6.7- to 12.2-GHz flux 
density ratio for these observations computed using the detection 
limit of the 12.2-GHz observations. 

a sample of water masers. Our observations show that of the 
four 6.7-GHz methanol masers in the LMC, two have an as- 
sociated water maser [IRAS and N160a/MC76), 
a detection rate of 50 per cent. Including the six new water 
masers we have detected, there are a total of thirteen water 
masers associated with star formation in the LMC. Hence 
the detection rate of 6.7-GHz methanol masers towards wa- 
ter masers is 15 per cent for the LMC, significantly lower 
than in the Galaxy and again consistent with the observed 
under-abundance of methanol masers in the LMC. 

4.2 Masers as star formation tracers 

We have searched the Spitzer SAGE Infrared Array Cam- 
era (IRAC) point source catalogue and archive for sources 
within 5 arcseconds of the maser positions. Four of the maser 
sources were associated (separation less than 5 arcsec) with 
sources in the SAGE IRAC catalogue. These are listed in 
Table (2] Somewhat surprisingly n one of these sou rces were 
identified as YSO candidates by IWhitnev et al. I dloos) m 
their catalogue of candidate young stellar objects. We have 
examined the Multiband Imaging Photometer (MIPS) 24- 
(Um images from the SAGE program and find that all the 
maser sources are associated either with a point source 
(listed in Table [2]), or are projected against more extended 
mid-infrared emission. 



iGruendl fc Chul l|2009l ) have independently searched for 
YSO candidates in the LMC. Their analysis, likewise, used 
the SAGE data, but also included a large amount of pub- 
lished and archival data from other wavelengths. We have 
compared the positions of the interstellar masers in the LMC 
with their catalogue and find that 11 of the 16 locations 
are within 2 arcseconds of an object they have classified as 
a "definite" YSO. The sources without an associated YSO 
candidate all lie in regions where there is a complex back- 
ground, or a bright source nearby (Robert Gruendl, private 
communication ) . It i s stri king that the methods used by 
IWhitnev et al] (|2008l ') and iGruendl fc ciiul l|2009l ) to iden- 
tify YSOs in the LMC (with the bulk of the observational 
data in common ) produ ce such very different results. Clearly 
IGruendl fc c"h^ (|2009D have been much more successful in 
identifying the young, high-mass objects with which masers 
are associated. 

We have compared the infrared properties of the maser- 
associated sources to those of the whole s ample of definite 
YSOs identified by ICruendl fc Chul (|2009l ) (their Table 9). 
Our initial examination of infrared colour-colour and colour- 
magnitude diagrams suggested that the maser-associated 
YSOs tended to be brighter in all near and mid-infrared 
wavelength bands and generally have redder colours than the 
bulk of YSO candidates (Figure (Tj). This is consistent with 
the maser-associated sources being more luminous and more 
deeply embedded than the bulk of YSO candidates, and may 
be useful for identifying regions of parameter space associ- 
ated with the rare, and highly sought youngest high-mass 
stars. The discrimination of the maser-associated sources is 
most clearly seen in the colour-magnitude diagrams, which 
explains why similar investigations undertaken with samples 
of Galactic masers typically produce less conclusive results. 
For Galactic sources the broad range of distances to the 
sources adds additional scatter to the observed magnitudes, 
an effect which is difficult to correct for when the distances to 
most sources are poorly known. For the [8.0] versus [4.5]- [8.0] 
colour-magnitude diagram the masers-associated sources all 
satisfy the criteria [4.5]-[8.0] > 2.5 and [8.0] < 12.5-([4.5]- 
[8.0]). A total of 352 of the 855 YSOs meet these criteria. 
For the [3.6] versus [3.6]-[4.5] colour-magnitude diagram the 
maser-associated sources (with one exception) satisfy the 
criteria [3.6]-[4.5] > 0.5 and [3.6] > 12.2-h([3.6]-[4.5]) and a 
total of 251 of the 855 YSOs in the LMC meet these crite- 
ria. The dashed lines in Fig.[7]indicate these criteria for each 
colour- magnitude diagram. Taking the union of the two sets 
results in a total of 195 YSOs satisfying both sets of criteria 
(22 percent of the sample as a whole). 

Recent statistical investigations of millimetre dust 
clumps with and without associated masers has enabled the 
development of models which predict the presence of a maser 
on the basis of the dust-clump properties (|Breen et al.ll2007l . 
[20091) . These models also allowed some inferences to be made 
about the critical physical properties which determine the 
presence of a maser. Motivated by the apparent differences 
between those YSOs with, and those without masers in the 
colour-magnitude diagrams we decided to undertake a sim- 
ilar inve stigation using the infr ared photometric data com- 
piled by IGruendl fc Chul (|2009l ') for the sample of YSOs in 
the LMC. We constructed boxplots for each of the nine wave- 
length bands available (J, H, K from the Two Micron All 
Sky Survey, 3.6, 4.5, 5.8 and 8.0 ^im from IRAC and 24 and 
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Table 2. The association of known interstellar masers in the LMC with infrared sources from SAGE. NOTE: In the N113/MC24 
field 051317.69-692225.0 is more than 2 arcseconds away from two of the masers, and we have given these separations merely as useful 
information. 
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Figure 7. Colour- magnitude diagrams comparing the sample of definite YSOs in the LMC from [Cruendl fc Chul Hooi) (black points), 
with the sub-sample of sources with an associated interstellar maser (red cross) and those sources within the maser fields without 
an associated maser (blue squares). The te n dency for the maser-associated sources to be associated with more luminous and redder 
sources can clearly be seen. iGruendl fc Chul l l200dh used a [4.5]-[8.0] versus [8.0] colour-magnitude diagram as the initial basis for their 
investigations. All the sources in their sample meet the criteria [4.5]-[8.0] > 2 and [8.0] < 14 - ([4.5]-[8.0]). The dashed lines in each plot 
indicate the region where most of the maser sources lie. For the left-hand plot this is [4.5]-[8.0] > 2.5 and [8.0] < 12.5-([4.5]-[8.0]) and 
for the right-hand plot this is [3.6]-[4.5] > 0.5 and [3.6] > 12.2-f ([3.6]-[4.5]). 



70 nm from MIPS) which show that in each of these wave- 
length bands the YSOs associated with masers tend to be 
more luminous tha n the sample as a whol e (Figure [8]). The 
photometric data in lGruendl fc Chul (|2009l ) is only complete 
for the 4.5- and 8.0-^m bands, as these were used in the ini- 
tial selection of YSO candidates. While most sources have 



measurements for all of the IRAC wavelength bands, this 
is not the case for the 2MASS and MIPS data, especially 
for the maser-associated YSOs. In order to objectively test 
whether the maser-associated YSOs are more luminous than 
the sample as a whole we performed a t-test for each wave- 
length band to determine if the two groups have a different 
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mean. The t-test returned a statistically significant result 
(> 95 per cent probability that the means of the two sam- 
ples are different) for each of the K-, 3.6-, 4.5-, 5.8-, 8.0- and 
24-/im bands which gave p-values of 0.03, 0.0003, 0.0004, 
0.0001, 0.0004 and 0.02 respectively. For the other wave- 
length bands the number of maser-associated YSO with a 
measurement was in all but one case, less than half the sam- 
ple (of 11 sources), which prevented meaningful statistical 
comparisons from being undertaken. 

We also investigated the infrared colours formed from 
the IRAC data. From the four wavelength bands six indepen- 
dent colours can be formed. The three colours which include 
data from the 3.6-/im band all have a statistically differ- 
ent mean (as measured by a t-test) for the maser-associated 
sources compared with the rest of the sample. The p- values 
for the difference in the mean for the [3.6]-[8.0], [3.6]-[5.8] 
and [3.6]- [4.5] colours are 0.03, 0.03 and 0.01 respectively. In 
each case the maser-associated sources have a larger mean 
colour (they are "redder"), than the rest of the sample . 

It would be desirable to fit a Binomial generalised linear 
model (GLM) to the infrared data to predict the properties 
of sources likely to be associated with a maser. However, in 
order to undertake such modeling the only sources which can 
be used are those where the presence or absence of a maser 
is well established, and which have data for all the parame- 
ters. For this reason, it was only practical to consider using 
the IRAC data for the GLM fitting as the data for the other 
wavelength bands are much less complete. Of the 8 55 definite 
YSO sources identified bv lCruendl fc Chul (|2009h 814 have 
observations for all four IRAC wavelength bands. Unfortu- 
nately, the small number of detected methanol masers and 
the targeted/incomplete nature of the water maser searches 
meant that for neither transition was GLM modeling prac- 
tical due to the small sample sizes. 

Fig.Quses blue squares to represent those YSO within a 
maser field which do not have an associated maser. It is clear 
that the majority of these sources have similar infrared prop- 
erties to the maser associated sources. This is in contrast to 
our findings when we compared the maser-associated YSOs 
to the rest of the sample. When we perform t-tests compar- 
ing those YSOs which fall within a maser field with the rest 
of the sample we find a statistically significant difference be- 
tween the means of the two samples for all the infrared bands 
with the exception of the MIPS 70 /im band (where there is 
insufficient data in the YSO near a maser sample). In addi- 
tion to this, the p- values from the f-tests are in all cases lower 
than when comparing the masers associated YSOs with the 
rest of the maser sample. This is likely to be because the 
presence of a strong water maser signposts a young high- 
mass star formation region and since they form in clusters 
there are likely to be other similar sources nearby. This is 
in contrast with most of the YSO sources, which differ sig- 
nificantly in their infrared properties from those associated 
with masers. They are likely to either be from more-evolved 
regions or associated with low- or intermediate-mass star 
formation. 

In summary, while we have established that there is a 
statistically significant difference between the infrared prop- 
erties of maser-associated YSOs compared to the whole sam- 
ple, the information we have on the presence or absence of 
masers is insufficient to enable us to undertake predictive 
statistical modeling at the present time. 



4.3 SED modeling 

Colour-colour and colour-magnitude diagrams are useful for 
identifying regions of parameter space which favour cer- 
tain types of objects. They are always subject to a de- 
gree of confusion however, due to overlaps in the colour 
ra nges of different types of objects (as nicely demonstrated 
by iGruendl fc C"h^l2009l ). In addition, they are often incom- 
plete (sources need measurements in each of three or four 
wavelength bands to be included), or they do not utilize all 
the available information (for sources where there are mea- 
surements in other wavebands). Fitting the spectral energy 
distribution (SED) of a source has the advantage that it can 
utilize all the available information for a p articular source. 
We ha ve used the online fitting SED tool of lRobitaille et al] 
(|2007l ) to investigate the SED's of those YSOs associated 
with masers and also the properties of the other LMC YSOs. 
The SED's are fitted to a grid of two-dimensional radiative 
transfer models of YSOs consisting of a central star sur- 
roun ded by an accreting disk and envelope (|Robitaille et al.l 
l2006f ). The limitations inherent in the assumptions, imple- 
mentation and results of the model have been extensively 
discu ssed in the literature (see for example lRobitaili3l2008l : 
I Grave fc Kumar 2009 ) and we will not repeat them here. 

In order to fi t each of the 855 definite YSO sources 
of IGruendl fc Chul l|2009l ). we used an automated script to 
populate the online fitting form with the infrared data and 
download the best-fit model parameters for each source. We 
have u sed the photometric measurements of lCruendl fc Chul 
(|2009h for each of the bands where they are available. When 
there was no reported measurement for a particular wave- 
band we have assumed upper limits of 17.6, 16.9 and 16.1 
magnitudes for the J-, H- an d Ks-band deep 2MASS ob- 
servations towards the LMC l|Cutri et al.ll2003 ) and upper 
limits of 19.3 and 16.1 magnitudes for the IRAC 3.6- and 5.8- 
/im bands (the observations are compl ete in the oth er IRAC 
bands due to the selection criteria) (jSewild [20081 ). Where 
no measurement was available for the MIPS 24- and IQ-^m 
bands we chose not to use the upper limit values in the fit, 
as in many cases the sources are in regions of extended or 
saturated emission at these wavelengths, rather tha n regions 
with no emission. Following IRobitaille et akl (|2007l ) we have 
assumed aperture sizes of 3 arcseconds for the 2MASS obser- 
vations, 5 arcseconds for IRAC observations, 10 arcseconds 
for MIPS 24 /im and 20 arcseconds for MIPS 70 /im. We 
allowed the distance to vary between 48 and 52 kpc and the 
visual extinction (along the line of sight) to vary between 
and 2 magnitudes (the mean value f or the interstellar exti nc- 
tion towards the LMC measured bv llmara fc Blit3 (|2007l 'l is 
0.3 magnitudes). Fig.[9]shows some examples of the resulting 
SED fits for maser and non-maser associated YSOs. 

There are between seven and nine infrared intensities 
or upper limits available for each source and in most cases 
the chi-squared for the resulting SED fits are reasonable. 
However, because of the large number of sources modeled 
we have made no attempt to remove sources where this was 
not the case. Figure [TO] shows the distribution of some of 
the SED model fit parameters for the YSOs associated with 
masers compared to the other YSOs in the lCruendl fc Chul 
(|2009l ) sample. For most of the parameters there is a large 
range present in both categories of YSO and in many cases 
the distributions are highly skewed (e.g. the envelope mass 



12 Ellingsen et al. 



1.25 [im 



1.65 i^m 



2.2 )im 






Other 



Maser 



Other 



Maser 



Other 



Maser 



3.6 )im 



4.5 )im 



5.8 )im 



— I — 
Other 



1 

Maser 



: X 



— I — 

Other 



1 

Maser 



— I 

Other 



1 

Maser 



8.0 )im 



24 |im 



70 |am 





CM 
I 



O - 
IM - 
^ - 



Other 



Maser 



Other 



Maser 



1 

Other 



1 

Maser 



Figure 8. Boxplots comparing tlie distribution of infrared intensity (magnitude) for various wavelength bands for sources known to be 
associated with masers and for the rest of the sample ("Other"). The y-axis for each box plot gives the intensity in magnitudes. The solid 
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mean of the maser and other distributions for all wavelength bands except 1.25 /im (J), 1.65 /im (H) and 70 fim. 



shown in Fig. I10|) . For skewed distributions it is more ap- 
propriate to use the Mann- Whitney test to determine if the 
two samples have a different median, rather than the t-test 
(which compares the means). The Mann- Whitney test re- 
turns statistically significant results (p- value < 0.05) for five 
of the model parameters, the mass of the central star (0.001), 
the outer radius of the envelope (0.05), the ambient density 



(0.01), the inclination to the line of sight (0.01) and the to- 
tal luminosity (0.0008). The inferences which can be drawn 
from these results is clear. The masers are associated with 
high-mass, luminous YSOs where the SED is dominated by a 
large envelope and are located in regions where the ambient 
density is higher than in the vicinity of most YSOs. 

The outer radius of the envelope and the inclination 
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Figure 9. The SED and fitted models for two maser associated sources (a & b), a source witii colours diff erent from mase r -assoc iated 
YSOs (c) and a YSO witfiin a maser field (d). The circles represent measurements taken from Table 9 of iGruendl &: Chul | |2009| ). the 
triangles are upper limits for 2MASS 6x and SAGE observations of the LMC. The best fit model is the solid black line and the grey 
lines show other acceptable models (chi-squared values within 3 of the best fit model). The dashed line shows the spectrum of the stellar 
photosphere as it would appear in the absence of circumstellar dust, but including interstellar extinction. The fitted mass of the central 
object and bolometric luminosity for each of these sources are (a) M=35.8M0 ; L = 1.26 X IO^Lq (b) M=28.9M0 ; L = 1.16 X lO^L© 
(c) M=8.4M0 ; L = 3.45 X IO^Lq (d) M=18.4Mq ; L = 9.34 X IO^Lq. 



to the line of sight the results appear to be influenced by 
details of the modeling process. The maximum outer ra- 
dius of the envelope used by the model is 10^ AU, and that 
is the fitted value for 9 of the 11 maser-associated YSOs 
(compared to 360 of the other 844 sources). Essentially this 
result is consistent with the masers being associated with 
higher mass objects, which would be expected to arise from 
the largest clumps/envelopes of gas and dust. Similarly for 
the inclination angle the maximum value available in the 
models is 87.13°, and that is the fitted value for 7 of the 
11 maser-associated YSOs (compared to 252 of the other 
844 sources). Naively, it seems unlikely that geometry of the 
star forming region should play a significant role in whether 
we see a luminous water maser associated with a particular 
source and such a result is likely to be an artifact of the 
modeling process. However, water masers are typically asso- 
ciated with outflows and an outflow-axis near perpendicular 
to the line of sight will produce the maximum line-of-sight 
velocity coher ent path length and hen ce the most luminous 
water masers. iRobitaiUe etal] l|2006h suggest that the in- 
clination angle is one of the most important parameters in 



determining the SED for young objects, as it controls the 
contribution from the inner, hotter regions of the envelope 
due to orientation of the bipolar cavity. The preference for 
the maser-associated YSOs to best fit models with high in- 
clination angles then means that there is little or no contri- 
bution from the inner regions in these SEDs. An alternative 
interpretation of these results would be that additional ob- 
scuring material is present for the highest mass YSOs in 
comparison to lower-mass objects. 

We have compared the SED model parameters of the 
maser-associated YSOs with other YSOs in the maser 
fields. We applied the Mann- Whitney test for all of the 
model parameters to compare these two samples and the 
only statistically significant result (p-value 0.03) was for the 
inclination to the line of sight. The fact that the inclination 
angle again shows a statistically significant difference war- 
rants further future investigation to determine whether it 
is an artifact of the modeling process, or if it has a physical 
basis. Since the SED model fits are based on the same 
infrared data for which we found no significant difference 
between the maser-associated YSOs and those in the maser 
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Figure 10. Boxplots comparing the distribution of some of the 
fitted parameters for the SED models for sources know to be asso- 
ciated with masers and for the rest of the sample ( "Other" ) . The 
six parameters shown here are the five for which a Mann- Whitney 
test shows a statistically significant difference in the median of the 
two distributions, plus the distribution of the envelope mass. See 
the caption for Fig. [S] for a description of how the boxplots can 
be interpreted. 



fields, it it not surprising that we see a similar lack of 
difference in the SED models. This means that, as for the 
infrared properties it is not possible to undertake GLM 
modeling to try and statistically predict which of the other 
YSOs in the LMC may have an associated luminous water 
maser. The best we can do is to look for other YSOs for 
which the fitted models are very similar to those associated 
with masers. Of the list of 844 YSOs in the LMC which 
are not known to have an associated water maser there 
are a total of 17 for which each of the mass of the central 
star, total luminosity, ambient density, outer radius of the 
envelope and inclination to the line of sight are all greater 
than or equal to the median of the maser-associated YSO 
sample. These YSOs are 044927.29-691206.3, 045358.57- 
691106.6, 045426.06-691102.3, 051004.00-685438.3, 
051328.26-692241.8, 051922.98-693940.5, 052208.52- 

675922.1, 052257.55-680414.1, 053205.89-662530.9, 

053603.03- 673234.2, 53831.62-690214.6, 053912.67-692941.4, 

053937.04- 694536.7, 053937.53-694609.8, 053945.94- 

693839.2, 054309.58-675006.6 & 054630.37-693458.6. 
Three of these (051004.00-685438.3, 051328.26-692241.8, 



053945.94-693839.2) are within one of the maser fields and 
so are known not to have an associated maser, however, 
the remaining 14 should be primary targets in any future 
searches for water masers in the LMC. 



5 CONCLUSION 

A sensitive search for 12.2-GHz methanol masers towards 
a sample of eight fields containing all currently known 
methanol, OH and water masers in the Large Magellanic 
Cloud has resulted in the detection of one source. This is 
the first detection of a 12.2-GHz methanol maser outside 
our Galaxy. Our sample size is too small and our obser- 
vations have insufficient sensitivity to determine if there is 
any difference in the properties of the different maser transi- 
tions in the LMC compared to the Galaxy; however they are 
consistent with the previously observed underabundance of 
methanol masers in the LMC. 

We have also observed the water masers associated with 
known high-mass star formation regions in the LMC and 
discovered six new sites. Five of these are in regions where 
emission had previously been detected, while the sixth is a 
weak new water maser detected towards the strongest 6.7- 
GHz methanol maser in the LMC IRAS 05011-6815. We find 
most of the star-formation masers in the LMC for which ac- 
curate positions are available are associated with sources 
identified as YSO on the basis of their infrared properties. 
The YSOs associated with masers tend to be more lumi- 
nous and redder than the sample as a whole. We suggest 
that this is likely due to the masers being associated with 
young high-mass star formation regions. SED modeling of 
the YSOs shows that the masers are preferentially associ- 
ated with sources with a high central mass, total luminosity 
and ambient density. This is consistent with expectations 
from Galactic maser studies. 
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